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r"| ' 1. Introduction 

' ' Understanding the physics of event shapes and jet rates (collectively jet 

observables) requires a deep knowledge of QCD dynamics. Their distribu- 
r~| ! tions (defined as the fraction of events S(w) with the observable's value V 

less than v) explore a wide range of physical scales, from the more inclusive 
. region v = 0{1), well described by fixed order perturbative (PT) calcula- 

I tions, to the extreme exclusive region t; — > where the language of quarks 

■ and gluons is no longer applicable. Describing the intermediate region re- 

quires a resummation of logarithmic enhanced contributions that appear at 
all orders in the PT expansion V. Resummation can then be seen as the 
link between PT and NP physics. Given that interest, a number of event 
shapes have been introduced, and for most of them there exist both next-to- 
leading order (NLO) calculations as well as resummed predictions at next 
to-leading logarithmic (NLL) accuracy. 

Unfortunately, while it is straightforward to compute a NLO distribution 
by just interfacing a fixed order Monte Carlo program with the observable 
definition in the form of a computer routine, for what concerns resummation 
a separate analytic calculation has to be performed for each observable. 

Together with Gavin Salam and Giulia Zanderighi we started then a 
project whose goal is the construction of a computer code that takes as 



* Talk presented at the XXXIII International Symposium on Multiparticle Dynamics, 
September 5-11, 2003, Krakow, Poland. 



(1) 



2 



BANFI'ISMD PRINTED ON FEBRUARY 1, 2008 



input from the user the definition of an observable and returns its resummed 
distribution at NLL accuracy. The starting point of such a programme is 
the classification of LL and NLL effects. 



Consider a Born system consisting of n hard partons {legs) {pi}i=i,...,n 
and an observable V that vanish in the Born limit, i.e. = 0. 

For most jet observables LL contributions can be addressed by comput- 
ing S(t') when a single soft gluon collinear to the hard legs is emitted and 
exponentiating the result. 

For any of these exponentiating quantities, NL logarithms can have dif- 
ferent origins: 

• emission of gluons soft at large angles or hard and collinear. Such 
contributions can be resummed by exponentiating the single emission 
result for T,{v); 

• multiple emission effects, which originate when, given a set of emitted 
partons {h}, all V{ki) < v but V{{ki}) > v (or viceversa), so that 
one has to take into account how all emissions contribute to build up 
the observable's value; 

• non-global logarithms ^2 , that arise whenever V is sensitive to emis- 
sions only in a part of the phase space. Unfortunately such NLL 
contributions are known only in the large Nc limit, thus reducing the 
accuracy of resummed predictions for non-global observables. 

Given the above classification we restrict ourselves to exponentiating 
global observables. Under the additional hypothesis that, after a soft emis- 
sion collinear to leg i, the observable's dependence on the emitted transverse 
momentum kt, rapidity r] and azimuthal angle (j) (all with respect to i) obeys 
the following parametrisation 



2. Automated resummation 
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at NLL accuracy T,{v) is given by the master formula [3] 






Here R{v) embodies all LL contributions and the NLL contributions that 
can be taken into account by exponentiating the single emission result. This 
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function, computed analytically, depends parametrically on a^, 6^, di as well 
as on the azimuthal average {higi{(j))). 

The NLL function J-{R') represents the multiple emission correction 
factor. Its general expression is |4j 

ni^) = (exp { -i?' In .yf/V'"^W^ U I) ' 

where the average is taken over all configuration of soft and collinear gluons 
such that, for fixed R' , the probability density of emissions along leg I is 
constant (the actual value of the constant depends on R' and on the colour 
factors of the hard legs, see 0). 

Although the above procedure is rather complicated, we were able to 
embody it a computer code, CAESAR, that completely automates the re- 
summation of any exponentiating global jet-observable j^]- 

3. Event-shapes in hadronic dijet production 

Let us show how the program works in the specific case of an event shape 
in hadronic dijet production. Given a unit vector n orthogonal to the beam 
axis we define the global transverse thrust: 

rt,, = 1 - max '^^ f , (4) 

where pu are the final state transverse momenta with respect to the beam. 

The program chooses a reference Born event with two incoming and two 
outgoing hard partons and generates soft and collinear emissions. It first 
checks that a resummation in the (2 + 2)-jet limit is feasible and that the 
observable is global. It then verifies that the parametrisation holds and 
determines for each leg ai,b£,d£. It recognises also whether ge{(j)) = \ sin(/)|'^'^, 
with Q integer, otherwise tabulates the azimuthal dependence. The program 
then should determine whether the observable exponentiates and compute 
the function J-{R') via a Monte Carlo procedure jlj. In this case it recognises 
that Tt^g is additive, that is V{ki^ . . . , kn) = ^(^i) + • • • + V{kn), which 
allows us to skip all the subsequent steps, since for additive observables 
exponentiation holds and J='{R') = e''^'^^' /r{l + R') [T]. 

These results are automatically tabulated (see fig.^), and used as inputs 
for the master formula © to compute the NLL resummed Tt^g distribution. 
In fig. n we see the differential distribution D{Tt^g) = d'E(Tt^g)/dlnTt^g at 
the Tevatron run II c.o.m. energy ^/s = 1.96TeV corresponding to selected 
dijets with Et > 50GeV and \ri\ < 1. We use the CTEQ6M pdf set and 
set factorisation and renormalisation scales at the partonic c.o.m. energy. 
The plot shows a clean separation among the various partonic channels, 
information that can be exploited for fits of the parton densities. 
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Fig. 1. A sample output of CAESAR: the results of the observable analysis (left) 
and the resummed distribution at NLL accuracy (right) for the global transverse 
thrust. The analysis has been performed with a reference Born configuration with 
the two outgoing partons back-to-back at an angle cos 6* — 0.2 with respect to 
the beam. The final result does not depend on the chosen configuration, as the 
program can check automatically. 



4. Conclusions and outlook 

The main feature of our approach is that the output curves are pure NLL 
functions, without contamination with spurious subleading contributions, so 
that they can be straightforwardly matched with fixed order results, and any 
hadronisation model can be applied. At present our effort is concentrated 
in releasing the first version of the program. The next step will be including 
in the code an automated matching with NLO calculations (6j, which will 
open the way to a vast amount of phenomenological analyses. 
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